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(N ■ ABSTRACT 

k>" ' We use a series of cosmological N-body simulations for a fiat ACDM cosmology to 

investigate the structural properties of dark matter haloes, at redshift zero, in the 
^ ■ mass range 3 x 10 9 h~ 1 M e < M vir < 3 x 10 13 /i -1 M Q . These properties include the 

iJ I concentration parameter, c, the spin parameter, A, and the mean axis ratio, q. For the 

concentration-mass relation we find c oc M^°' n in agreement with the model proposed 
\q ' by Bullock et al. (2001), but inconsistent with the alternative model of Eke, Navarro 

, & Steinmetz (2001). The normalization of the concentration-mass relation, however, 

is 15 percent lower than suggested by Bullock et al. (2001). The results for A and 
q are in good agreement with previous studies, when extrapolated to the lower halo 
masses probed here, while c and A are anti-correlated, in that high-spin haloes have, 
Q ■ on average, lower concentrations. In an attempt to remove unrelaxed haloes from the 

Jjj | sample, we compute for each halo the offset parameter, x g, defined as the distance 

c/3 . between the most bound particle and the center of mass, in units of the virial radius. 

Removing haloes with large x s increases the mean concentration by ~ 10 percent, 
lowers the mean spin parameter by ~ 15 percent, and removes the most prolate haloes. 
In addition, it largely removes the anti-correlation between c and A, though not en- 
tirely. We also investigate the relation between halo properties and their large-scale 
environment density For low mass haloes we find that more concentrated haloes live 
in denser environments than their less concentrated counterparts of the same mass, 
consistent with recent correlation function analyses. Note, however, that the trend is 
weak compared to the scatter. For the halo spin parameters we find no environment 
dependence, while there is a weak indication that the most spherical haloes reside in 
slightly denser environments. Finally, using a simple model for disk galaxy formation 
we show that haloes that host low surface brightness galaxies are expected to be hosted 
by a biased sub-set of haloes. Not only do these haloes have spin parameters that are 
larger than average, they also have concentration parameters that are ~ 15 percent 
lower than the average at a given halo mass. We discuss the implications of all these 
findings for the claimed disagreement between halo concentrations inferred from LSB 
rotation curves, and those expected for a ACDM cosmology. 

Key words: galaxies: haloes - cosmology:theory, dark matter, gravitation - methods: 
numerical, N-body simulation 
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1 INTRODUCTION 

The theory of cold dark matter (CDM) provides a success- 
ful framework for understanding structure formation in the 
universe. Within this paradigm dark matter collapses first 
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into small haloes which merge to form progressively larger 
haloes over time. Galaxies are thought to form out of gas 
which cools and collapses to the centers of these dark mat- 
ter haloes (White & Rees 1978). 

In the standard picture of disk galaxy formation the 
structural and dynamical properties of disk galaxies are ex- 
pected to be strongly related to the properties of the dark 



2 A. V. Maccid et al. 



matter haloes in which they are embedded. In particular 
the characteristic sizes and rotation velocities of disk galax- 
ies are determined (to first order) by the spin parameter, 
concentration parameter, and size of the host dark matter 
halo (e.g. Mo, Mao & White 1998, hereafter MMW). Conse- 
quently, the detailed rotation curve shapes of disk galaxies 
can, in principle, be used to constrain the structural prop- 
erties of their dark matter haloes. This is especially true for 
low surface brightness (LSB) galaxies, which are believed to 
be dark matter dominated even at small radii. A steadily 
increasing data base of observed LSB rotation curves has 
resulted in a heated debate as to whether the slopes of the 
inner density profiles of dark matter dominated disk galaxies 
are consistent with the cuspy profiles found in N-body simu- 
lations, or similarly, whether the inferred concentrations are 
as high as predicted (see Swaters et al. 2003 and references 
therein). 

Unfortunately, determining cusp slopes and/or concen- 
tration parameters from mass modeling rotation curves is 
non- unique, even for dark matter dominated galaxies (e.g., 
Dutton et al. 2005). In particular, determining c requires 
knowledge of the virial radius, which is hard to constrain 
using data that only covers the inner ~ 10% of the halo. 
As an alternative measure of the central density of a halo, 
Alam, Bullock & Weinberg (2002) introduced a dimension- 
less quantity that does not require knowledge of the halo 
virial radius, and demonstrated convincingly that the ob- 
served rotation curves of LSB galaxies imply halo concen- 
trations that are systematically lower than predicted for a 
flat ACDM cosmology with a matter density Qm = 0.3 and 
a scale-invariant Harrison-Zeldovich power-spectrum with 
normalization ag = 1.0. 

Further observational support for a lower normaliza- 
tion of the c — M v j r relation comes from the zero point of 
the rotation velocity-luminosity relation, also known as the 
Tully-Fisher relation (Tully & Fisher 1977) of disk galaxies 
(van den Bosch 2000). Detailed disk formation models have 
clearly demonstrated that the high concentrations of CDM 
haloes cause an overprediction of the rotation velocities at 
a fixed disk luminosity, at least for a 'standard' ACDM cos- 
mology with SIm = 0.3, Oa = 0.7, and as = 0.9 (e.g., Dutton 
et al. 2007; Gnedin et al. 2006). 

Although these discrepancies may indicate a genuine 
problem for the CDM paradigm, there are a number of al- 
ternative explanations: First of all, as shown by various au- 
thors (e.g. Swaters et al. 2003; Rhee et al. 2004; Spekkens, 
Giovanelli, & Haynes 2005) the observed rotation curves 
could be hampered by a variety of observational biases, such 
as beam smearing, slit offsets and inclination effects, all of 
which tend to underestimate the circular velocity in the cen- 
tral regions. 

Secondly, the dark matter distribution could have been 
modified by astrophysical processes such as bars (e.g. Holley- 
Bockelmann, Weinberg, & Katz 2005) or dynamical fric- 
tion (e.g. Mo & Mao 2004; Tonini et al. 2006). These pro- 
cesses tend to lower the concentration of the dark matter 
halo, bringing it in better agreement with the observations. 
On the other hand, adiabatic contraction (e.g., Blumenthal 
et al. 1986) thought to be associated with the formation of 
disk galaxies, actually tends to increase the halo concentra- 
tion, and it remains to be seen whether the above mentioned 
processes are strong enough to undo this contraction and 



still cause a relative expansion of the inner halo (see Dutton 
et al. 2007 for a detailed discussion) . 

A third option is that the data-model comparison has 
been made for the wrong cosmology. In particular, a re- 
duction in the power of cosmological density fluctuations 
on small scales causes a significant reduction of the pre- 
dicted halo concentrations (e.g. Eke, Navarro & Steinmetz 
2001; Zentner & Bullock 2002; Alam et al. 2002; van den 
Bosch, Mo, & Yang 2003). Most data-model comparisons 
have been based on a flat ACDM cosmology with Qm = 0.3 
and as = 0.9. However, recently the third year data re- 
lease from the WMAP mission has advocated a model with 
flu =± 0.25 and cr 8 ~ 0.75 (Spergel et al. 2006). This rel- 
atively small change in cosmological parameters causes a 
significant reduction of the predicted halo concentration pa- 
rameters, bringing them in much better agreement with the 
data (e.g., Yang, Mo, & van den Bosch 2003; van den Bosch, 
Mo & Yang 2003). 

Another potentially important cause for the discrep- 
ancy are systematic errors in the actual model predictions. 
Both Bullock et al. (2001a; hereafter B01) and Eke, Navarro 
& Steinmetz (2001; hereafter ENS) presented analytical 
models that allow one to compute the mean halo concen- 
tration for given halo mass, redshift and cosmology. Unfor- 
tunately, at redshift zero the predictions of these models are 
divergent below ~ W 11 Mq, with the ENS model predict- 
ing halo concentrations that are significantly lower. This is 
of particular importance for LSB (and dwarf) galaxies with 
V < 100 km s" 1 , which are thought to typically reside in 
haloes with masses below this value. Both B01 and ENS cal- 
ibrated their models against numerical simulations. Those 
of B01 probed the mass range between 3 x 10 Mq and 
6x 10 13 /i _1 M Q , while those of ENS probed an even narrower 
range from 3 x 10 11 /i- 1 M e to 3 x 10 12 /i _1 M© (albeit with 
higher resolution). What is needed to discriminate between 
these models are simulations that resolve a large population 
of low mass haloes, which is one of the main objectives of 
this paper. 

Another important issue that we wish to address in this 
paper is the possibility that the LSB disk galaxies that have 
been used to constrain halo concentrations reside in a biased 
sub-set of haloes (see discussion in Wechsler et al. 2006). Nu- 
merical simulations have shown that there is a significant 
scatter in both halo concentration, c, and halo spin parame- 
ter, A, at a given halo mass (e.g. Bullock et al. 2001a,b). Thus 
if disk galaxies form in a biased subset of haloes, this could 
lead to an apparent discrepancy between theory and obser- 
vation. In fact, there are a number of potential causes for 
such a bias. First of all, disk galaxies are expected to prefer- 
entially form in haloes that have not experienced any recent 
major merger. There is evidence that such a subset of haloes 
has higher mean c, lower mean A, and lower scatter in both 
c and A (Wechsler et al. 2002; D'Onghia & Burkert 2004). 
Clearly, this would worsen the disagreement between model 
and data. On the other hand, it has also been suggested that 
LSB galaxies preferentially reside in haloes with relatively 
low concentrations. First of all, since disks are thought to 
be in centrifugal equilibrium, less concentrated haloes will 
harbor less concentrated (i.e., lower surface brightness) disk 
galaxies (e.g., MMW, B01). In addition, using numerical 
simulations Bailin et al. (2005) found that haloes with higher 
spin parameters have, on average, lower concentration pa- 
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rameters. Since LSB galaxies are thought to be those with 
high spin parameters such a correlation would imply that 
LSB galaxies reside in haloes with relatively low concentra- 
tions. If confirmed this could offer an alternative explanation 
as to why (some) LSB galaxies have lower concentrations 
than predicted. Note, however, that previous studies (BOf , 
Navarro, Frenk & White 1997, hereafter NFW), have found 
no correlation between spin parameter and concentration. 

Another potential bias for disk galaxy formation could 
arise if there is a correlation between environment (defined 
as the large scale matter density) and c or A. In particular, 
Harker et al. (2006) found evidence that low mass haloes in 
dense environments assemble earlier than haloes of the same 
mass in under-dense environments. Note however, that for 
the lowest density environments this trend reverses so that 
formation redshifts actually increase with decreasing den- 
sity. Since haloes that assemble later are less concentrated 
(Wechsler et al. 2002) , one thus may expect a similar correla- 
tion between halo concentration and environmental density. 
If dwarf and LSB galaxies preferentially form in under-dense 
regions this could also help explain the lower than expected 
halo concentrations of these galaxies. 

In this paper we study galaxy size dark matter haloes 
from a set of cosmological N-body simulations with the fol- 
lowing goals: (i) to test the predictions of B01 and ENS 
regarding the halo concentrations of low mass haloes (down 
to ~ 3 x 10 9 /i _1 Mq), (ii) to determine the scatter in con- 
centration, spin parameter, and halo shape at a given mass, 
(iii) to determine whether there is a correlation between the 
spin and concentration parameters, and (iv) to determine 
whether c, A and halo shape depend on the density of the 
environment in which the halo is located. Our paper is or- 
ganized as follows. Section 2 describes our set of N-body 
simulations. In Section 3 we discuss how halo concentra- 
tion, halo spin parameter and halo shape depend on halo 
mass. Sections 4 investigates whether c and A are correlated, 
while Section 5 focuses on the environment dependence of 
halo properties. In Section 6 we use simple models for disk 
formation to investigate whether one expects LSB galaxies 
to reside in haloes with a biased concentration parameter. 
Finally, we summarize our results in Section 7. 



2 N-BODY SIMULATIONS 

In order to explore as wide a range of virial masses as pos- 
sible, we run simulations of 4 different box sizes, listed in 
Table 1. For comparison we also show the parameters of the 
Bullock et al. (2001a,b), Bailin et al. (2005) and millennium 
run (Springel et al. 2005) simulations. The B01 simulation 
has similar size and mass resolution as our 64 a ,b boxes, while 
our smallest box simulation has a mass resolution that is ~ 5 
times higher than that of Bailin et al. (2005). In order to test 
for cosmic variance, and to increase the size of our sample 
we ran two simulations for each of the three smallest box 
sizes. 

All simulations have been performed with PKDGRAV, 
a tree code written by Joachim Stadel and Thomas Quinn 
(Stadel 2001). The code uses spline kernel softening, for 
which the forces become completely Newtonian at 2 soften- 
ing lengths. Individual time steps for each particle are cho- 
sen proportional to the square root of the softening length, e, 



Table 1. N-body Simulation Parameters 



Name Box size N particle mass force soft. 



[ft- 1 Mpc] [/i _1 M©] [h-^pc] 

14a, 14b 14.2 250 3 1.4 X 10 7 0.43 

28a, 28b 28.4 250 3 1.1 x 10 8 0.85 

64a, 64b 63.9 300 3 7.2 x 10 8 1.92 

128 127.8 300 3 5.8 x 10 9 3.83 



Bullock 60 256 3 1.1 x 10 9 1.8 

Bailin 50 512 3 7.8 x 10 7 5.0 

Millennium 500 2160 3 8.6 X 10 8 5.0 



over the acceleration, a: Aij = r\^j e/a,i. Throughout, we set 
r\ = 0.2, and we keep the value of the softening length con- 
stant in comoving coordinates during each run. The physical 
values of e at z = are listed in Table 1. Forces are computed 
using terms up to hexadecapole order and a node-opening 
angle 8 which we change from 0.55 initially to 0.7 at z = 2. 
This allows a higher force accuracy when the mass distri- 
bution is nearly smooth and the relative force errors can be 
large. 

We adopt a flat ACDM cosmology with parameters 
from the first year WMAP results (Spergel et al. 2003): 
matter density J7m = 0.268, baryon density fij, = 0.044, 
Hubble constant h = H /(W0 kms" 1 Mpc -1 ) = 0.71, and 
a scale-invariant, Harrison-Zeldovich power-spectrum with 
normalization ag = O.fl The initial conditions are gener- 
ated with the GRAFIC2 package (Bertschinger 2001), which 
also computes the transfer function as described in Ma & 
Bertschinger (1995). The starting redshifts Zi are set to the 
time when the standard deviation of the smallest density 
fluctuations resolved within the simulation box reaches 0.2 
(the smallest scale resolved within the initial conditions is 
defined as twice the intra-particle distance, while the maxi- 
mum scale is set by the Box size). 

In all of our numerical simulations, haloes are iden- 
tified using a SO (Spherical Overdensity) algorithm. As a 
first step, candidate haloes are located using the standard 
friends-of-friends (FOF) algorithm, with a linking length 
fen -1 / 3 , with n the mean particle density and b a free pa- 
rameter which we set to 0.2. We only keep FOF haloes with 
at least iV m t n = 200 particles, which we subject to the fol- 
lowing two operations: (i) we find the point, C, where the 
gravitational potential due to the group of particles is mini- 
mum, and (ii) we determine the radius 7? v ir, centered on C, 
inside of which the density contrast is A v j r . For our adopted 
cosmology A v j r ~ 98 (using the fitting function of Mainini 
et al. 2003). Using all particles in the corresponding sphere 
we iterate the above procedure until we converge onto a sta- 
ble particle set. The set is discarded if, at some stage, the 
sphere contains less than A^min particles. If a particle is a 
potential member of two haloes it is assigned to the more 
massive one. For each stable particle set we obtain the virial 
radius, i? v ir, the number of particles within the virial radius, 

1 The recent analysis of the three year WMAP data (Spergel 
et al. 2006) suggests lower values for Cm, <?8 and the spectral 
index. In a forthcoming paper (Maccio et al. , in preparation) 
we investigate the effects of these new cosmological parameter on 
our results. The main change regards a lower normalization of the 
concentration, as expected from the B01 and ENS models. 
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N v i r , and the virial mass, Af v i r . Above a mass threshold of 
7Vvir = 250 particles there are ~ 2750, 3750, 7450 and 4500 
haloes in the simulations of box size 14.2, 28.4, 63.9, 127.8 
h~ 1 Mpc respectively (these numbers refer to the two ver- 
sions of each box size combined together). 

In Fig. [T] we report the comparison of the halo mass 
functions of all our simulations with the analytical mass 
function of Sheth & Tormen (2002). Since the Sheth & Tor- 
men mass function has been tuned to reproduce the mass 
function of FOF haloes (with b — 0.2), we use the same FOF 
masses here. For the remainder of this paper, however, we 
consistently will use the spherical overdensity masses, M v i r , 
described above. Note that the FOF mass functions agree 
well with the Sheth & Tormen mass function over the full 
five orders of magnitude in halo mass probed by our simu- 
lations: all data points are consistent with the model within 
one sigma (error bars show the Poisson noise in each bin due 
to the finite number of haloes). Moreover all the simulations 
made with different box sizes agree with each other in the 
mass ranges where they overlap. 

2.1 Halo parameters 

For each halo we determine a set of parameters as described 
below. All of these parameters are derived from the SO 
haloes (i.e. from the particle sets defined by the SO criteria), 
rather than from the FOF particle sets. 



2.1.1 concentration parameter 

N-body simulations have shown that the spherically aver- 
aged density profiles of DM haloes can be well described by 
a two parameter NFW profile: 



(1) 



P(r) 5c 

Pcrit (r/r s )(l + r/r B ) 2 ' 

where p cr it is the critical density of the universe, 8 C is the 
characteristic overdensity of the halo, and r s is the ra- 
dius where the logarithmic slope of the halo density profile 
dlnp/dlnr = —2 (NFW). A more useful parametrization is 
in terms of the virial mass, M v j r , and concentration param- 
eter, c = .Rvir/Vs. The virial mass and radius are related by 
M v i r = A v i r /9 cr it(47r/3)-R vir , where A v i r is the density con- 
trast of the halo. 

To compute the concentration of a halo we first deter- 
mine its density profile. The halo center is defined as the 
location of the most bound halo particle, and we compute 
the density (pi) in 50 spherical shells, spaced equally in log 
radius. Errors on the density are computed from the Pois- 
son noise due to the finite number of particles in each mass 
shell. The resulting density profile is fit with a NFW profile 
(Eq. [1} , which provides a good fit to most haloes over the 
range of radii we are interested in. Note that, in this paper, 
we are not concerned with the inner asymptotic slope of the 
density profile. During the fitting procedure we treat both 
r s and S c as free parameters. Their values, and associated 
uncertainties, are obtained via a x 2 minimization procedure 
using the Levenberg & Marquart method. We define the 
r.m.s. of the fit as: 



JV 



where p m is the fitted NFW density distribution. We do 
not use the \ 2 value of the best-fit since this increases with 
TVvir. This occurs because higher resolution haloes have bet- 
ter resolved substructure and smaller Poisson errors on the 
density, thus making the fit worse. Finally, we compute the 
concentration of the halo, c, using the virial radius obtained 
from the SO algorithm, and we define the error on log c as 
(<r rs /r s )/ ln(10), where er rs is the fitting uncertainty on r s . 

We checked our concentration fit pipeline against the 
one suggested by B01. As a test we used both the procedures 
to compute c in all our cubes. No systematic offset arises in 
the concentration vs. mass relation due to the different halo 
definition and fitting procedure. 



2.1.2 spin parameter 

The spin parameter is a dimensionless measure of the 
amount of rotation of a dark matter halo. The standard 
definition of the spin parameter, due to Peebles (1969), is 
given by 

7.I/5. 1 1/2 

'J vir -L^vir 



A = 



cm:; 



5/2 



(3) 



where M v i r , Jvir and _E v ir are the mass, total angular mo- 
mentum and energy of the halo, respectively. Due to difficul- 
ties with accurately measuring J5 v ir, Bullock et al. (2001b) 
introduced a modified spin parameter: 



A' 



</vir 



(4) 



(2) 



with Vvir the circular velocity at the virial radius. For a sin- 
gular isothermal sphere these two definitions are equivalent. 
For a pure NFW halo, however, they are related according 
to A = A'/(c) 1/2 with /(c) = §c[(l + cf - 1 - 2(1 + c) ln(l + 
e)]/[c-(l + c)ln(l + e)] 2 (MMW). In what follows, we define 
A^ = A'/(c) as the "corrected" spin parameter. 

In order to avoid potential problems and inaccuracies 
with the measurement of E v i T , we adopt the A' definition 
for the halo spin parameter, unless specifically stated oth- 
erwise. The advantage of A' over X' c is that the latter can 
introduce artificial correlations between halo concentration 
and halo spin parameter, since an error in c translates 
into an error in /(c) and hence A c . We define the uncer- 
tainty in log A' as (aj /J) / ln( 10), were we use that aj/J — 
+ 0.04/A' 2 ) ~ 0.2/\'yf{N) (Bullock et al. 2001b). 
Note that this implies that the errors on A' are largest for 
haloes with a low spin parameter and with few particles. 



2.1.3 shape parameter 

Determining the shape of a three-dimensional distribution 
of particles is a non-trivial task (i.e. Jing & Suto 2002) . Fol- 
lowing Allgood et al. (2006) we determine the shape of our 
haloes starting from the inertia tensor. As a first step the 
inertia tensor of the halo is computed using all the particles 
within the virial radius; in this way we obtain a 3 x 3 ma- 
trix. Then the inertia tensor is diagonalized and the particle 
distribution is rotated according to the eigen vectors. In this 
new frame (in which the moment of inertia tensor is diago- 
nal) the ratios 02/01 and 03/01 (01,02,03 being the major, 
intermediate and minor axis, respectively) are given by: 



Concentration, Spin and Shape of Dark Matter Haloes 5 



02 
«1 



'E^M^i 

E 



03 
01 



E'< 



(5) 



Next we again compute the inertia tensor, but this time 
only using the particles inside the ellipsoid defined by 01, a 2 
and 03. When deforming the ellipsoidal volume of the halo, 
we keep the longest axis (01) equal to the original radius 
of the spherical volume (cf. Allgood et al. 2006) . We iterate 
this procedure until we converge to a stable set of the axis 
ratios. 

Since dark matter haloes tend to be prolate, a useful 
parameter that describes the shape of the halo is q = (0,2 + 
03) /2oi, with the limiting cases being a sphere (q = 1) and 
a needle (q = 0). 



2.1.4 offset parameter 

The last quantity that we compute for each halo is the offset, 
Xoff, defined as the distance between the most bound particle 
(used as the center for the density profile) and the center of 
mass of the halo, in units of the virial radius. This offset is a 
measure for the extent to which the halo is relaxed: relaxed 
haloes in equilibrium will have a smooth, radially symmet- 
ric density distribution, and thus an offset that is virtually 
equal to zero. Unrelaxed haloes, such as those that have only 
recently experienced a major merger, are likely to reveal a 
strongly asymmetric mass distribution, and thus a relatively 
large x b . Although some unrelaxed haloes may have a small 
x Q g, the advantage of this parameter over, for example, the 
actual virial ratio, 2T/V, as a function of radius (Maccio, 
Murante & Bonometto 2003; Shaw et al. 2006), is that the 
former is trivial to evaluate. 

Fig. [2] shows histograms of, and correlations between 
iVvir, Prms, x Q g, and q. The rms of the density profile fit 
decreases with iVvin as expected, while x a s and q axe uncor- 
rected with iVvir (especially for N V - 1T > 10 3 ). The solid lines 
in the x ff-N v it plot show the ratios of the softening length 
to the virial radius. This shows that the offset parameter is 
not effected by resolution effects. In an ideal halo simulation 
with a large number of particles x r is expected to decrease 
with the decrease of the halo mass. We do not see this trend 
in the XoB-Nvir relation because in our simulations the value 
of x s is mostly dominated by numerical effects at the low 
mass tail. Moreover since we used different simulations with 
different mass resolution there is not a one to one relation 
between mass and AT V ir- Note also that x D g is uncorrelated 
with p rms , but that there is a strong correlation between x s 
and q so that more prolate haloes tend to have larger offsets. 
We discuss the implications of this correlation in £|3 . 1 1 

The distributions of p rms and log(a; g) are approxi- 
mately normal with means of 0.4 and —1.4, respectively. The 
distribution of q, on the other hand, is strongly skewed. Note 
also that the ^-distributions are slightly different for differ- 
ent box-sizes. This is because there is a correlation between 
halo mass and halo shape, such that more massive haloes 
are less spherical (see Section ETT1 below) . Since larger boxes 
contain more massive haloes, this results in lower mean axis 
ratios. 



3 MASS DEPENDENCE OF SPIN, 
CONCENTRATION AND SHAPE 

Fig. [3] shows the relations of concentration vs. mass, and 
spin parameter vs. mass for all haloes with more than 250 
particles within the virial radius. The right panels show the 
data from the simulations with the four box sizes clearly 
visible. The left panels show the mean (solid dots) and twice 
the standard deviation of c and A' (error bars) in bins equally 
spaced in logAf v i r (plotted at the mean M V [ T in each bin). 
The mean c and A' are computed taking account of both 
the estimated measurement errors and the intrinsic scatter, 
using 

V = . Wi 



,2,2- 
("hit + fi . 



(6) 



Here t/i denotes either c or A' of the ith halo, Wi is the weight 
(haloes with larger uncertainties receive less weight), and a t 
is the measurement error on c or A'. The intrinsic scatter 
(Tint is given by 



fint = 



£» [(yi - yY 



(7) 



Since these two equations are coupled the computation of y 
and (Ji n t requires an iterative procedure. We start by com- 
puting y and (Ji nt with Wi = 1. Next we iterate until a stable 
solution for y and a lnt is found. This procedure converges 
rapidly, typically in 3-4 iterations. 

The solid red lines in Fig.[3]show weighted (using errors 
on c and A') least squares fits of c and A' on Af v i r . The 
c — A/ V i r relation is well fitted by a single power-law over 
four orders of magnitude in mass 3 x 10 9 ft -1 Mq < M v i r < 
3 x lO 13 h _1 M 0) with 

logc= 1.020[±0.015] - 0.109[±0.005](log M viT - 12). (8) 

Note that M vlr is in units of h~ Mq . The numbers in square 
brackets give the scatter in the corresponding value between 
the seven different simulations. The total scatter about this 
mean relation is <ri nc = 0.40±0.03 and the intrinsic scatter is 
fine — 0.33 ±0.03, where again the uncertainty corresponds 
to the box-to-box scatter. These are in excellent agreement 
with the total and intrinsic scatter found by B01 which are 
0.41 and 0.32, respectively (see also Wechsler et al. 2002). 
The slope of the A' — M v j r relation is consistent with zero, in 
agreement with previous studies (e.g. Lemson & Kauffmann 
1999, hereafter LK99; Shaw et al. 2006). If we take A' to be 
independent of M v i r , we find a median A' = 0.034 ± 0.001 
with an intrinsic scatter of a in v = 0.55 ± 0.01, which are in 
excellent agreement with Bullock et al. (2001b) (median of 
A' = 0.035±0.005 and a lnX , = 0.50±0.03) and other studies 
(e.g., van den Bosch et al. 2002; Avila- Reese et al. 20050- 

The dashed and dotted lines show the mean c(M v i r ) for 
the B01 and ENS models, respectively. In addition, for both 
models we also show the upper and lower 2a intrinsic scatter 
bounds, where we adopt ai nc = 0.32 (Wechsler et al. 2002). 



2 Haloes with lower A' have, on average, larger errors, which re- 
sults in a tail of haloes with low A'. As long as the larger mea- 
surement uncertainties of these haloes are properly taken into 
account, this tail does not effect the mean and scatter of the A' 
distribution. 
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The B01 model has two free parameters: K, which deter- 
mines the normalization of the relation, and F, which effects 
the slope. The ENS model has just one free parameter C a 
which controls the normalization; the slope is fully deter- 
mined by the model. If we adopt F = 0.01, as advocated 
by B01, we find that the slope of the B01 model is in excel- 
lent agreement with our simulations over the full range of 
masses covered. Consequently, our data is inconsistent with 
the significantly shallower slope of the ENS model. Note that 
ENS compared their model to a very small sample of relaxed 
haloes (~ 15), albeit with high resolution (iV V ir > 30000), 
over a small mass range M vil ~ 2 x 10 11 — 4 x 10 12 h" 1 M . 
Over this mass range the ENS model is in reasonable agree- 
ment with our simulation data. 

In terms of the normalization our data is best fit with 
K — 3.4 ± 0.1 (for F = 0.01), where the error reflects 
the effect of cosmic variance as determined from the vari- 
ous independent simulation boxes. This is 15 percent lower 
than the K = 4.0 originally advocated by B01, but consis- 
tent with Zhao et al. (2003) and Kuhlen et al. (2005), who 
found K — 3.5 (both for F = 0.01). In their paper Kuhlen 
et al. argued that the cause for their lower normalization 
might be due to the A^-body code used for their simulations 
(GADGET-1). However, we have used an independent code, 
PKDGRAV, and obtain the same result. We therefore sus- 
pect that the cause for this discrepancy resides somewhere 
else. Indeed, as it turns out, B01 used a slightly different 
transfer function for setting up the initial conditions of their 
numerical simulation than for computing their model. If they 
correct this, they obtain a best fit K — 3.75 (James Bullock 
& Andrew Zentner, private communication). If we take our 
(admittedly crude) estimate for the cosmic variance at face 
value, this still implies that we predict a normalization that 
is significantly lower (at the 3<r level). 

As pointed out by B01, a better match to the slope 
of the c — M v i r relation for haloes more massive than ~ 
10 13 /i -1 Mq (at the expense of worsening the agreement for 
haloes with M < lO 11 /! -1 !^©) can be obtained by using 
F — 0.001, in which case we find that K — 2.6. This is 
again approximately 15 percent smaller than the K — 3.0 
advocated by B01 for this value of F. 

3.1 The Impact of Unrelaxed Haloes 

Our halo finder (and halo finders in general) do not distin- 
guish between relaxed and unrelaxed haloes. There are two 
reasons why we might want to remove unrelaxed halos. First, 
and most importantly, unrelaxed haloes often have poorly 
defined centers, which makes the determination of a radial 
density profile, and hence of the concentration parameter, an 
ill-defined problem. Secondly for applications to disk galaxy 
formation, haloes that are not in dynamical equilibrium are 
unlikely to host disk galaxies, and even less likely to host the 
more dynamically fragile LSB galaxies. In this case, the halo 
parameters inferred from (LSB) disk rotation curves need to 
be compared to those of the subset of relaxed haloes only. 

One could imagine using p rms (the r.m.s. of the NFW 
fit to the density profile) to decide whether a halo is relaxed 
or not. However, while it is true that p rms is typically high 
for unrelaxed haloes, haloes with relatively few particles also 
have a high p rma (due to Poisson noise) even when they are 
relaxed. This is evident from Fig. [2] which shows that p rm s 



and iVvk are strongly correlated. Furthermore not all unre- 
laxed haloes have a high p rms . We found several examples of 
haloes with low p rms which are clearly unrelaxed. This is due 
to the smoothing effects of spherical averaging when com- 
puting the density profile. However, these haloes are often 
characterized by a large x g (the offset between the most 
bound particle and the center of mass) . In what follows we 
therefore use both p rms and x Q g to judge whether a halo is 
relaxed or not. 

Fig. [4] shows the residuals of the c — M v i r and A' — M v i r 
relations for haloes with iVvir > 250, against iVvir, Prms, x Q g, 
and q. The filled circles and error bars show the mean and 2a 
scatter of points in equally spaced bins. The smaller error 
bars, sometimes barely visible, show the Poisson error on 
the mean (cr/yN). Both the c and A' residuals show no 
trend with iV v j r , even down to our limit of 250 particles. 
This indicates that numerical resolution is not affecting our 
results. 

Interestingly we find that haloes with the lowest p rms 
tend to have larger c, smaller scatter in c, and lower A'. 
A similar result for halo c was found by Jing (2000), who 
used the maximum deviation of the density profile from the 
NFW fit, which is similar to our p rms parameter. In Table 2 
we show the effect of the RMS parameter on the mean and 
scatter of halo c and A'. Here < logCi2 > is the zero point 
of the c — M v i r relation, measured at M V1I = 10 12 /i _1 Mq. 
To allow for comparison with Jing (2000) we only use haloes 
with iVvir > 10000. We find that the highest p rm s haloes have 
a mean c which is roughly a factor of 2 lower, and a scatter 
in c which is roughly a factor of 2 higher, than the lowest 
prms haloes, in agreement with Jing (2000). Additionally we 
also find a factor of 2 difference between the average A' of 
the highest and lowest p rm s haloes. Haloes with the highest 
Prms have the highest mean offset parameter, < x g >, and 
the lowest mean halo shape parameter, < q >. This suggests 
that these haloes are the most unrelaxed. 

The trends between c and A' with x g are shown in the 
third column of Fig. [3] These show that for haloes with small 
x g the residuals are uncorrelated with x s- However, these 
haloes have concentrations that are higher and spin param- 
eters that are lower than the overall average. For x a g > 0.04 
(log x g > —1.4) there is a clear systematic trend that haloes 
with larger x a have lower c and higher A'. The same trends 
are seen for q, where more prolate haloes (of a given mass) 
have lower c and higher A'. This basically reflects the corre- 
lation between x g and q (see Fig. [2}. 

These trends are consistent with unrelaxed haloes being 
the systems that experienced a recent major merger: (i) the 
center of the halo is poorly defined, which results in a large 
x B an d an artificially shallow (low concentration) density 
profile, (ii) the halo shape is more prolate, and (iii) the spin 
parameter is higher due to the orbital angular momentum 
'transferred' to the system due to the merger (e.g., Vitvit- 
ska et al. 2002; Mailer, Dekel & Somerville 2002). Ideally 
one would test the correspondence between merger histories 
and x g using the actual merger trees extracted from the 
simulation. Unfortunately, we did not store sufficient out- 
puts to be able to do so. We intend to address these issues 
in a future paper, based on a new set of simulations. 

We now split our sample into 4 sub-samples according 
to prms and x a ff, with dividers of p rm s = 0.4 and x Q g — 
0.04, which correspond to the mean of the distributions of 
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Table 2. Effect of RMS parameter on c and A' distributions for haloes with iV vrf > 10000. 
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prms 


N haloes 


< log C12 > 


°"lnc 


< log A' > 




< Prms > 


< a=ofi > 


< q > 


> 0.25 


88 


0.770 


0.392 


-1.285 


0.528 


0.36 


0.100 


0.52 


0.15-0.25 


183 


0.986 


0.281 


-1.436 


0.596 


0.19 


0.038 


0.66 


0.10-0.15 


226 


1.111 


0.204 


-1.552 


0.513 


0.12 


0.022 


0.73 


0.00-0.10 


117 


1.138 


0.167 


-1.625 


0.516 


0.09 


0.018 


0.74 


> 0.0 


614 


1.039 


0.353 


-1.493 


0.597 


0.17 


0.031 


0.68 



p rms and x ff, respectively (see Fig. [2]|. We refer to the 4- 
subsamples as: 

• GOOD (p rms < 0.4 & x oS < 0.04); 

• BAD (p mis < 0.4 & x oS > 0.04); 

• UGLY (p rms > 0.4 & x oS > 0.04); 

• NOISY (p rms > 0.4 & XoS < 0.04). 

Fig. [5] shows the c — M v i r and A' — M v j r relations for the 
haloes in the GOOD sub-sample with iVvir > 250. As with 
the full sample, this c — M v i r relation is well fitted with a 
single power-law given by 

logc = 1.071 [±0.027] - 0.098[±0.009](log M vir - 12). (9) 

This relation has a slope that is ~ 10 percent shallower and 
a zero point that is ~ 10 percent higher than for the full 
sample (eq. [8]). The total scatter about this mean relation 
is <Ti llc = 0.30 ± 0.03 and the intrinsic scatter is <ri nc = 
0.26 ±0.03, about 25 percent lower than for the full sample. 
The B01 model again accurately fits the c — M v i r relation, 
but with K = 3.7 ± 0.15 (for F = 0.01). For F = 0.001 the 
best-fit value for K is 2.9. 

The slope of the c(M) relation for low mass haloes is still 
consistent with that predicted by the B01 model but steeper 
than the prediction of the ENS model. Thus, the fact that 
ENS only compared their model to a small sample of well 
relaxed haloes can not explain the discrepancy between their 
results and those of B01. As already eluded to above, the 
main reason the ENS model was found to be consistent with 
their own simulations is that the haloes in their simulation 
only covered a small range in halo masses, over which the 
difference in slope with respect to the B01 model is difficult 
to infer. 

The slope of the A' — M v i r relation is again consis- 
tent with zero. However, the median is ~ 15 percent lower 
(A' = 0.030 ± 0.003) and the intrinsic scatter is reduced 
by ~ 5 percent (cr lnA ; = 0.52 ± 0.01). These differences in 
c(Mvir) between the full set of haloes and our GOOD sub- 
sample are almost identical to those obtained by Wechsler 
et al. (2002) between all haloes and haloes without a ma- 
jor merger since z = 2. This reinforces the notion that our 
GOOD sub-sample consists mostly of haloes which have not 
experienced a recent major merger. 

Fig. [S] shows the dependence of the halo shape parame- 
ter q on halo mass. We find that more massive haloes are on 
average more flattened (more prolate), in qualitative agree- 
ment with previous studies (Jing & Suto 2002; Kasun & 
Evrard 2005; Bailin & Steinmetz 2005; Allgood et al. 2006). 
As shown in Allgood et al. (2006), the halo shape is fairly 
tightly related to the halo assembly time, such that haloes 
that assemble later are less spherical (and less relaxed). To 
first order this explains the decrease of q with increasing 



halo mass, as well as the relation between q and x g (see 
Fig. [2]| . Fig. [6] also shows that haloes with high x Q g (BAD 
and UGLY sub-samples) have a lower median q, as well as 
a much more pronounced tail to low q. Note also that there 
are very few highly prolate haloes (q < 0.5) with low x a s 
(NOISY and GOOD sub-samples). A potentially important 
implication of this is that (LSB) disk galaxies, which are too 
fragile to survive in unrelaxed haloes, are unlikely to reside 
in strongly flattened haloes. 



4 CORRELATION BETWEEN SPIN AND 
CONCENTRATION 

In their study of the halo angular momentum distributions, 
Bullock et al. (2001b) noticed that haloes with high A' have 
concentration parameters that are slightly lower than av- 
erage. Although such a correlation may not be totally un- 
expected, since both the spin parameter and concentration 
parameter depend on the mass accretion history (Wechsler 
et al. 2002; Vitvitska et al. 2002), Bullock et al. argued that 
the c-A' anti-correlation is a mere 'artifact' from the fact that 
(i) c and A are uncorrelated (see also NFW), and (ii) A and 
A' are related via /(c) (see Section [2.1.2[) . However, Bailin 
et al. (2005) used the X' c definition for the spin parameter 
(which should be equal to A) and claimed a significant anti- 
correlation between halo concentration and spin parameter. 
As discussed in Bailin et al. (2005), such an anti-correlation 
may have important implications for the interpretation of 
the rotation curves of LSB disk galaxies (see discussion in 
Section [T]). Using our large suite of simulations, we therefore 
re-investigate this issue. 

Bailin et al. (2005) focused on haloes with 1 x 
10 11 h~ 1 M e < Mvir < 2 x lO 12 ^- 1 !/© and iVpop > 1290. 
In Fig. [7] we plot the halo concentration as a function of the 
three different definitions of the spin parameter; A', X' c and A. 
To facilitate the comparison with Bailin et al. we only select 
haloes in our simulations with 1 x 10 11 }i~ 1 Mq < M v i r < 
2 x 10 12 h~ 1 M e and iV vir > 1000. All three plots reveal a 
weak but significant correlation in that the lowest concen- 
tration haloes have relatively low spin parameters. Contrary 
to NFW and Bullock et al. (2001b), but in agreement with 
Bailin et al. (2005) , we therefore argue that c and A are cor- 
related. Note that resolution is not an issue here, since we 
obtain the same result in each of our different simulations. 

Due to the mass dependence of c, a more illustrative way 
to look for a correlation between halo concentration and spin 
parameter is to plot the residuals (at constant M v i r ) of the 
c — M v i r and A' — M v i r relations. This is shown in Fig. [S]for 
all haloes with N v i r > 250. The lower left panel shows the 
residuals for the full sample, which reveals a clear trend that 
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haloes with high A' have lower than average c. We now split 
the haloes into the four sub-samples defined above, and plot 
their residuals (with respect to the mean c(M) and mean A' 
for the GOOD sub-sample). This shows that the correlation 
between c and A' is, at least partially, due to the inclusion of 
haloes with high x g (i.e., haloes that are unrelaxed), inde- 
pendent of whether the halo has a high or low p rms : the cor- 
relation is clearly more pronounced for the UGLY and BAD 
sub-samples. However, the GOOD and NOISY sub-samples 
still reveal a small trend that haloes with larger A' offsets 
have a lower c offset. We have repeated this analysis us- 
ing A^ and A, and the small correlation between residuals in 
the GOOD and NOISY sub-samples remains. We therefore 
conclude that there indeed is a small intrinsic correlation 
between c and A' at a fixed M v j r . However, when exclud- 
ing the unrelaxed haloes, the amplitude of this correlation 
is very weak compared to the scatter in both parameters. 
In Section [B] below we investigate to what extent this small 
correlation may effect (LSB) disk galaxies. 



5 ENVIRONMENT DEPENDENCE 

We now investigate whether the concentration, spin param- 
eter and shape of dark matter haloes are correlated with 
the large scale environment in which they are located. This 
is interesting because disk galaxies, and LSBs in particular, 
are preferentially found in regions of intermediate to low 
density (Mo et al. 1994; Rosenbaum & Bomans 2004). If 
haloes in low density environments have different structural 
properties than haloes of the same mass in over-dense en- 
vironments, this would therefore imply that the haloes of 
(LSB) disk galaxies are not a fair representation of the av- 
erage halo population. 

Using a set of numerical simulations, for different cos- 
mologies, LK99 studied the environment dependence of dark 
matter haloes. The only halo parameter that was found to 
be correlated with environment is halo mass (denser envi- 
ronments contain more massive haloes). Halo concentration, 
spin parameter, shape and assembly redshif10 were all found 
to be uncorrelated with halo environment. 

However, using exactly the same simulations as LK99, 
Sheth & Tormen (2004) found clear evidence that halo as- 
sembly times correlate with environment based on a marked 
correlation function analysis. A similar result was obtained 
by Gao, Springel & White (2005), who used the millennium 
simulation (Springel et al. 2005) to demonstrate that the 
low mass haloes that assemble early are more strongly clus- 
tered than haloes of the same mass that assemble later. For 
massive haloes, however, the dependence on the assembly 
time was found to diminish. This has since been confirmed 
by a number of studies (Harker et al. 2006; Zhu et al. 2006; 
Wetzel et al. 2007; Jing, Suto & Mo 2007). In addition, it 
has been found that the clustering of haloes also depends 
on other halo properties, such as halo concentration, halo 
spin parameter, subhalo properties, and the time since the 
last major merger (Wechsler et al. 2006; Zhu et al. 2006; 
Wetzel et al. 2007; Gao & White 2007; but see also Percival 

3 The assembly time of a halo of mass Mq is defined as the time 
at which the most massive progenitor has a mass equal to Mq /2. 



et al. 2003), which is not too surprising given that all these 
properties correlate with halo formation time (e.g., Wech- 
sler et al. 2002; van den Bosch, Tormen & Giocoli 2005). 
All these results seem to overrule the finding by LK99, and 
seem to suggest that dark matter properties depend rather 
sensitively on their large scale environment. 

Using our large sample of objects, we investigate 
whether the concentration, spin parameter and shape of 
haloes are correlated with the large scale environment in 
which they are located. Rather than using (marked) cor- 
relation functions, we follow LK99 and correlate the halo 
properties with the overdensity, 1 + S r , measured in spheres 
of radius r(with r = 1, 2, 4, 8 h _1 Mpc) centered on each halo: 

1 + Sr= <P(r)> = 3M«r) 1 

p n 47rr3 pcritfiM 

The results for all haloes are shown in Fig. [§] The top 
row shows the relation between virial mass and overden- 
sity. The vertical dashed line shows the overdensity at the 
virial radius. The solid diagonal line indicates the mass scale 
M = |7r(l+5 r )p cr i t nMr 3 . Thus all haloes with R v i r less than 
the filter radius should lie below this line, as is the case. Note 
that we do not include sub- haloes in our analysis, so the only 
haloes that can have very high densities log (I + 8) > 2 must 
be haloes with virial radii close to the filter radius, as is 
the case. We see that more massive haloes tend to live in 
more overdense regions (cf. LK99). As the filter radius is 
increased the mean overdensity tends towards zero, and the 
scatter in overdensity is reduced, as expected. Note that we 
do not compute the overdensity on 8/i _1 Mpc scales in the 
simulations with the smallest box size. 

The second, third and fourth rows of Fig. [5] show the 
residuals (at fixed M v i r ) of the c-M v i r , A'-M v j r and q-M v i r 
relations, respectively, all as function of overdensity. The 
roughly horizontal lines with errorbars indicate the mean 
residual (plus its errors) as function of overdensity, while 
the dashed lines outline the ±lcr scatter. Note that, within 
the errors on the mean, there is no significant indication that 
the residuals are larger or smaller in regions that are over- 
or under-dense. In other words, these results are in excellent 
agreement with those of LK99 and seem to suggest that halo 
properties are independent of their large scale environment. 
Note that our results cover a much wider dynamic range 
in halo mass, and are based on higher-resolution simula- 
tions than in the case of LK99. How can this be reconciled 
with the findings based on the correlation function analyses 
described above? Some insight is provided by the roughly 
vertical lines with errorbars (in the second, third and fourth 
rows of Fig. [9j , which indicate the average overdensity (plus 
its error) for haloes in a given residual bin. The correspond- 
ing dotted lines outline the ±I<r scatter. These show that 
haloes with the largest concentration, largest spin parame- 
ter, and/or that are most spherical (all with respect to the 
average for their mass) are located in slightly denser envi- 
ronments. Since the correlation function reflects ensemble 
averages, this seems consistent with the findings that more 
concentrated haloes and/or those with a larger spin parame- 
ter are more strongly clustered. We emphasize, though, that 
the trends seen in Fig. [5] are (i) weaker when environment is 
measured over a larger volume, (ii) only reveal an environ- 
ment dependence at the extremes of the residual distribu- 
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tions, and (iii) largely vanish when we only focus on haloes 
in our GOOD sub-sample (not shown). 

As emphasized in Harker et al. (2006) , in order to reveal 
the correlation between environment and assembly time, it 
is important to only consider haloes in a relatively small 
mass range. Therefore, in Fig. QJJ] we plot the residuals of 
the c — M v i r relation versus overdensity in 2h _1 Mpc spheres 
for haloes in mass bins with a width of 0.5 dex. The vari- 
ous lines and errorbars have the same meaning as in Fig. [5] 
For haloes with M v i r ^ 10 12 hT 1 Mq there is a weak, mildly 
significant indication that haloes in denser environments 
have slightly higher concentrations (reflected by close-to- 
horizontal lines). More pronounced, however, is the trend 
that more concentrated haloes live on average in denser 
environments (reflected by close-to- vertical lines). This is 
consistent with the fact that more concentrated, low mass 
haloes are more strongly clustered (cf. Wechsler et al. 2006). 
ft is clear from Fig. [TO] that the environment dependence is 
weaker for more massive haloes; in fact for haloes in the 
mass range 12 < log [M/(/i _1 M Q )] < 13.5 no significant 
environment dependence is seen (this explains why the sig- 
nal is weaker in Fig. [9] where we add all the masses to- 
gether). Again this is in good agreement with the clustering 
results of Wechsler et al. (2006) and Gao & White (2007), 
who found that the dependence on halo concentration van- 
ishes for haloes with masses close to the typical collapsing 
mass, M*: for the cosmology assumed in our simulations 
M* — 6.7x 10 12 hT 1 Mq. Thus, our results are in good agree- 
ment with the various claims based on correlation function 
analyses. They also illustrate, though, that the trends are 
weak compared to the scatter, which explains why LK99 
did not notice any environment dependence. Only when one 
carefully estimates the average environment as function of 
halo property, which is what a correlation function measures, 
does the environment dependence reveal itself. 

Figure [TT1 shows the residuals, at fixed M v i r , of the A' — 
Mvir relation as function of log(l+<$2) for the same mass bins 
as in Fig. 1101 Unlike the concentration, the spin parameter 
reveals no significant environment dependence, in any mass 
bin. This is inconsistent with Gao & White (2007) who found 
that higher spin haloes are more strongly clustered than low 
spin haloes of the same mass. It is unclear why our analysis 
of the environmental densities does not recover this trend, 
while it does reproduce the trend for the halo concentrations. 

Finally, Fig. [12] shows the same as Figs [10] and [TTJ but 
for the shape parameter q. As for the spin parameter, there 
is no indication for any significant environment depencence. 
This seems at odds with Fig. [9] (fourth row) , which shows 
that the most spherical haloes (those with large positive 
Ag) live in denser environments. This apparent discrepancy 
simply reflects number statistics: only the haloes in the bin 
with Aq > 0.2 seem to reside in regions that are denser than 
average. When split in mass bins, however, there are too few 
haloes with Ag > 0.2 to reveal the signal. 



6 THE HOST HALOES OF LSB DISK 
GALAXIES 

We now investigate whether LSB disk galaxies are expected 
to reside in dark matter haloes that form a biased sub-set in 
terms of their concentration parameters. In the MMW model 



the central surface density of an exponential disk, Eo,d, is 
determined by A, c, M V [ T , and the galaxy mass fraction m ga i 
(defined as the ratio between disk mass and halo mass). A 
lower Eo,d will result from 

• a higher A at fixed c, M V [ T , and m ga i; 

• a lower M V [ T at fixed c, A, and ra ga i; 

• a lower c at fixed A, M v i r , and m ga i; 

• a lower m ga i at fixed c, A, and M v i r . 

To complicate matters, lower mass haloes have higher con- 
centrations, A and c are weakly anti-correlated (see Sec- 
tion 2J, and m ga i is expected to increase with M v i r due to 
various astrophysical feedback processes. 

To investigate the impact of all these relations on the 
surface brightness of disk galaxies, we construct MMW type 
models as described in Dutton et al. (2007). These models 
consist of an exponential disk, a Hernquist bulge and a NFW 
halo. The halo parameters are M v i r , c and A c (= A), which we 
take directly from the haloes of our GOOD sub-sample. An 
additional parameter is the galaxy mass fraction m ga i , which 
we fix to m ga i = 0.04 for simplicity. The bulge formation 
recipe is based on disk instability, and therefore only effects 
the highest surface brightness disks; the details of this bulge 
formation recipe are not important for this work. We assume 
that the halo is unaffected by the formation of the disk, i.e. 
we do not consider adiabatic contraction. As highlighted in 
Dutton et al. (2007), models with adiabatic contraction are 
unable to simultaneously match the zero-point of the Tully- 
Fisher relation and the galaxy luminosity function. 

The MMW formalism gives the galaxy mass, M ga i, 
baryonic disk scale length, _Rd, and central surface den- 
sity of the baryonic disk, £o,d- As described in Dutton 
et al. (2007), we split the disk into a stellar and a gaseous 
component assuming that all disk material with a surface 
density E(J?) > Ecrit(-R) nas been turned into stars. Here 
E cr it is the star formation threshold density, modeled as the 
critical surface density given by Toomre's stability criterion 
(Toomre 1964; Kennicutt 1989). The resulting stellar surface 
density profile is fitted with an exponential profile to obtain 
the scale length of the stellar disk, and the central sur- 
face density of the stellar disk, Eo,*. Note that in general 
R* < R d - 

In order to facilitate a comparison with observations we 
convert Eo.* into /zo,b, the central surface brightness of the 
stellar disk in the B-band, using the B-band stellar mass-to- 
light ratio, Tb- Using data from Dutton et al. (2007) and 
relations between mass-to-light ratios and color from Bell 
et al. (2003) we obtain 



logTs = 0.06 + 0.25 log 



Alt 



1O 1O M 



(11) 



where we have adopted a Kennicutt IMF and a Hubble con- 
stant h = 0.7. In principle this relation has a scatter of ~ 0.1 
dex, which we ignore for simplicity. 

Fig. 1131 shows correlations between c, Ac, Eo.d (central 
surface density of baryonic disk), Eo,» (central surface den- 
sity of stellar disk), and ho,b (central surface brightness 
in the B-band) and histograms of Eo,d, Eo,* and /io,s for 
GOOD haloes in the mass range lO lo /i _1 M0 < M vir < 
10 13 /i -1 Mq. The distributions of £o,d, Eo,* and ^io,b are 
approximately log-normal, reflecting the log-normal distri- 
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bution of \' c . The distribution of /xo,b has a peak value in 
agreement with the Freeman value of 21.65 magn. arcsec -2 . 

As expected, there is a strong correlation between Eo,d 
and X' c , in that haloes with larger spin parameters host lower 
surface density disks. The flattening of this relation at low 
X' c owes to our bulge formation recipe. At fixed A^, more 
massive haloes have higher £o,d, despite their (on average) 
smaller concentrations. The scatter in c at a fixed Al v i r re- 
sults in some overlap between the three mass samples, but 
the three mass ranges are clearly visible. Thus, the depen- 
dence of surface density on halo mass is at least as important 
as the dependence on halo concentration. The same trends 
are seen in the relation between surface density of the stel- 
lar disk and spin parameter. However the mass separation is 
no longer present in the relation between surface brightness 
and spin parameter. This is because, at a given £o,*, more 
massive haloes have higher stellar mass-to-light ratios, and 
hence lower surface brightness. 

The lower left panel of Fig. [13] plots the halo concen- 
tration versus the halo spin parameter. Although these pa- 
rameters seem uncorrelated, there is a weak anti-correlation 
between c and A^ , as discussed in Section [4] The fact that 
this correlation is not as pronounced as in Fig. [7] is due to 
the fact that here we only consider the GOOD sub-sample. 
The other three panels in the lower row show the correla- 
tion between halo concentration and disk surface density (or 
brightness) . For haloes of a given mass, there is a clear corre- 
lation in that more concentrated haloes host higher surface 
density (brightness) disks. This correlation has two origins: 
centrifugal equilibrium and the (weak) correlation between 
A^ and c. In what follows we investigate the relative impor- 
tance of both of these causes. 

Fig. [T5] shows the distribution of c for different ranges 
of M v i r and /xo,b • The three mass ranges roughly correspond 
to massive galaxies (150 < V V i T ^ 300 km s -1 ), interme- 
diate mass galaxies (70 < V v i T < 150 km s -1 ), and dwarf 
galaxies (30 < V v \ t < 70kms -1 ). At a fixed surface bright- 
ness, less massive haloes have higher c as expected. At a 
fixed halo mass, there is a clear trend that lower surface 
brightness disks reside in less concentrated haloes. If we de- 
fine LSB galaxies as those with a central surface brightness 
24 > Ho,b ^ 23 magn. arcsec -2 , we find that they live in a 
sub-set of haloes whose average concentration (at fixed halo 
mass) is ~ 15 percent lower than the overall average for 
that halo mass. We therefore conclude halo concentrations 
inferred from LSB rotation curves should not be compared 
to (c)m, but rather to /(c) a/, with / ~ 0.85 a bias cor- 
rection factor. A similar conclusion was obtained by Bailin 
et al. (2005), except that they found a bias correction factor 
of / ~ 0.70. As discussed in Section [4] this owes to the fact 
that they did not remove unrelaxed haloes from their sam- 
ple. Since we consider it unlikely that LSB galaxies reside in 
unrelaxed haloes, we belief our bias correction factor to be 
more realistic. 

Finally, in order to investigate the origin of this bias, 
we have run a control sample with the same distributions 
of M v i r , c and X' c as the simulation data, but with no corre- 
lation between c and A^. This reduces the correction factor 
to / ~ 0.95, and therefore shows that the main contribu- 
tion to / owes to the (very weak) anti-correlation between 
halo concentration and halo spin parameter. The remain- 
ing contribution simply reflects that centrifugally supported 



disks in less concentrated haloes will be less concentrated 
themselves. 



7 SUMMARY 

In this paper we have used a set of cosmological N-body 
simulations to study the scaling relations, at redshift zero, 
between the concentration parameter, c, spin parameter, A', 
shape parameter, q, and mass, M v i r , of a large sample of 
dark matter haloes. Due to the combined set of simula- 
tions, we were able to extend previous studies to a mass 
range at least an order of magnitude smaller, covering the 
full range of masses in which galaxies are expected to form: 
3 x lO 9 /! -1 M < M vir < 3 x lO^h" 1 M . 

For this mass range we find c oc M~ { °' , which is in 
agreement with the model of Bullock et al. (2001a), but in 
disagreement with the model of Eke, Navarro & Steinmetz 
(2001) which predicts a significantly shallower slope. The 
single free parameter of the ENS model only controls the 
normalization of the c — M v ir relation, so that their model 
cannot be tuned to fit the data. The ENS model has also 
been shown to be unable to match the slope of the c — M v i r 
relation for low mass haloes at z — 3 (Colm et al. 2004). For 
the Bullock et al. (2001a) model our data is well fitted with 
a model with F = 0.01 and K — 3.4 ±0.1 (where the error 
reflects our estimate of cosmic variance) . Note that this nor- 
malization is 15 percent lower than the K — 4.0 originally 
proposed by Bullock et al. (2001a), but it is in good agree- 
ment with Zhao et al. (2003) and Kuhlen et al. (2005), who 
found a best-fit normalization of K — 3.5. This discrepancy 
is at least partially due to an inconsistency with the use of 
transfer functions in the work of Bullock et al. If they use 
the same transfer functions to set up the initial conditions 
of their simulations and to compute the model predictions, 
they obtain K = 3.75. This is however significantly higher 
than can be accounted for by our estimate of cosmic vari- 
ance. We find an intrinsic scatter in c and A' at fixed M v i r 
of u\nc ~ 0.33 ± 0.03 and <Ji n \ — 0.55 ± 0.01, and a median 
spin parameter of A' = 0.034 ±0.001, all in good agreement 
with Bullock et al. (2001a,b). 

In an attempt to distinguish between relaxed and unre- 
laxed haloes we introduce a new and simple parameter: x Q g 
which is defined as the distance between the most bound 
particle and the center of mass, in units of the virial radius. 
The distribution of x fi is approximately log-normal with a 
median x g ~ 0.04. The full set of haloes shows strong cor- 
relations between x ff and the residuals of the c — M V [ T and 
A' — A/ V i r relations, such that haloes with larger x ft have a 
larger than average A' and a lower than average c. Removing 
haloes with large x a therefore results in a higher mean con- 
centration, a lower mean spin parameter, and in less scatter 
in both the c—M v [ T and A' — A/ V j r relations. The median spin 
parameter of 'relaxed' (GOOD) haloes is A' = 0.030 ± 0.003 
with an intrinsic scatter of <Ji n y = 0.52 ± 0.01. The average 
c(M v j r ) of 'relaxed' haloes is again in good agreement with 
the model of B01, but with F = 0.01 and K = 3.7, and with 
a reduced intrinsic scatter of (7\ nc = 0.26 ± 0.03. 

A better fit to the c— M v i r relation for high mass haloes 
(M vir > 10 13 /i _1 M ) can be obtained with F = 0.001 and 
K = 2.6 (for all haloes), and K = 2.9 (for GOOD haloes). 
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This is, however, at the expense of under predicting the con- 
centrations for the low mass haloes (M vir < 10 n /i _1 Mq). 

We also find a strong correlation between the mean axis 
ratio of the haloes, q, and x a s, such that more prolate haloes 
(i.e., those with lower q) have higher x g, on average. This 
suggests that the majority of haloes with small axis ratios 
are unrelaxed, and thus that fragile LSB galaxies are un- 
likely to reside in haloes that are strongly aspherical. This 
makes it less likely that the discrepancy between observed 
and predicted rotation curves is due to the fact that disks 
are strongly elliptical, as suggested by Hayashi et al. (2004). 

We have also investigated the environment dependence 
of the residuals in c and A'. This is interesting because disk 
galaxies, and LSBs in particular, are preferentially found 
in regions of intermediate to low density (Mo et al. 1994; 
Rosenbaum & Bomans 2004). Defining 'environment' by 
the matter overdensity within spheres of radii 1,2,4 and 
8h~ Mpc, we find at the low mass end (M < M*) that 
more concentrated haloes live in denser environments than 
their less concentrated counterparts of the same mass. This 
is consistent with the studies of Wechsler et al. (2006) and 
Gao & White (2007) who found a similar trend using cor- 
relation functions. However, as we have shown, this trend 
is weak compared to the scatter, which explains why Lem- 
son & Kauffmann (1999) did not notice any environment 
dependence in their simulation. Contrary to the halo con- 
centrations, the halo spin parameters reveal no environment 
dependence at fixed mass. This is at odds with the results of 
Gao & White (2007) who found that haloes with a large spin 
parameter are more strongly clustered than low spin haloes 
of the same mass. Finally we find a weak trend that the most 
spherical haloes reside in slightly denser environments. 

Finally, using a simple model for disk galaxy forma- 
tion, we investigated the properties of the (expected) host 
haloes of LSB disk galaxies (i.e., those with a central surface 
brightness 24 > /io,b > 23 magn.arcsec -2 ). In addition to 
having higher than average spin parameters, in agreement 
with numerous other studies (e.g., Dalcanton et al. 1997; 
Jimenez et al. 1998), we also find that the host haloes of LSB 
galaxies have concentrations that are biased low by about 
15 percent. This correlation between halo concentration and 
disk surface brightness (or density) owes largely to a (weak) 
anti-correlation between A' and c, and to the fact that cen- 
trifugal equilibrium commands that less concentrated haloes 
host less concentrated disks. The amplitude of this correla- 
tion is significantly smaller than what has been advocated 
by Bailin et al. (2005), but this owes to the fact that these 
authors did not remove unrelaxed haloes, which are unlikely 
to host fragile LSB galaxies, from their analysis. 

All these results have important implications for the in- 
terpretation of the halo concentrations inferred from LSB ro- 
tation curves. Numerous studies in the past have argued that 
these are too low compared to the predictions for a ACDM 
cosmology (e.g., Alam et al. 2002; McGaugh et al. 2003; de 
Blok et al. 2003). However, there are several reasons why 
we now believe that the model predictions where too high. 
First of all, virtually all previous predictions were made for 
a flat ACDM cosmology with Q m = 0.3 and erg = 0.9 (or 
erg = 1.0). However, if one adopts the cosmology favored by 
the three year data release of the WMAP mission (Spergel 
et al. 2006), one predicts concentrations that are about 25 
percent lower (Maccio et al. , in preparation). Compared to 



a ACDM model with Q m = 0.3 and ag — 1.0 the concentra- 
tions are 35 percent lower. Secondly, the B01 model (with 
F — 0.01 and K = 4.0) overpredicts the halo concentra- 
tions by ~ 15 percent. If we take into account that LSBs 
only reside in relaxed haloes, this is lowered to a ~ 8 per- 
cent effect. And finally, one needs to correct for the fact that 
LSB galaxies reside in a biased sub-set of haloes, which is 
another 15 percent effect. Combining all these effects, it is 
clear that the halo concentrations predicted were almost a 
factor two too large. This brings models and data into much 
better agreement. 
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Figure 1. Mass functions for the simulations. The point and color types correspond to the different runs: 14a (cyan triangles); 14b (red 
triangles); 28a (blue squares); 28b (yellow squares); 64a (magenta open circles); 64b (green open circles); 128 (black circles). The solid 
line is the Sheth & Tormen prediction for erg = 0.9. 
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Figure 2. Correlations between, 2V v ir> Prms, % tf > an d <?• The color coding in the density plots is according to the logarithm of the number 
of points in each cell. For the histograms the colors correspond to the different simulations as in Fig. [T] the thick black line shows the 
histograms of the combined samples. The black lines in the x a f{ — N v y plot shows the ratio of the softening length to the virial radius. 
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Figure 3. Relations between concentration, spin parameter and virial mass for haloes with N v - Ir > 250. The dashed lines give the mean 
and 2a scatter from Bullock et al. 2001a, b. Dotted lines give the ENS model prediction. The solid red lines give weighted power-law fits: 
y = zero + slope(log M v i r — 12), where y = log c or log A and M v i r is in units of /i"'Mq . The parameters of the fits are given in the lower 
left corner of each panel. The panels on the left show the data binned in mass. The filled circles give the error weighted mean, the small 
error bars gives the Poisson error on the mean while the larger error bars give the intrinsic 2-sigma scatter. The panels on the right show 
all the data points color coded according to the simulation: 14a (cyan); 14b (red); 28a (blue); 28b (yellow); 64a (magenta); 64b (green); 
128 (black). 
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Figure 4. Residuals of the c — M v i r (upper panels) and A' — M v i r (lower panels) relations as fitted in Fig. [3] against iV v i r , p r ms, x a g, 
and q for all haloes with N v „ > 250. The large error bars show twice the la intrinsic scatter, while the small error bar shows the possion 
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Figure 8. Correlations between the residuals of the c — M vlr and A' — M v - lr - relations in Fig. \5\ The lower left panel shows all haloes, 
while the remaining panels show the four sub-samples defined according to p rms and x a g as indicated. The error bars show the 2<r scatter 
in c for each AA' bin. 
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Figure 9. Environment dependence of M v jx, and of the residuals, at fixed Af v i r of the c — M v i T , A' — M v - lr - and q — M v i r relations. Here 
environment is defined by the matter overdensity within spheres of radii 1,2,4 and 8 h~ 1 Mpc. In the top row the diagonal line shows the 
minimum overdensity for a halo of a given mass, and the vertical dashed line shows the overdensity at the virial radius. In each panel, 
the solid curve shows the mean, the error bars indicate the Poisson errors on the mean, and the dashed curves show the ±lcr scatter in 
each mass bin. In the second, third and fourth rows the almost vertical lines indicate the average overdensity (plus its error) for haloes 
in a given residual bin. 



Concentration, Spin and Shape of Dark Matter Haloes 21 




0.5 


-0.5 





_) 










1 ' 





I I I I I I I I I I I I I I 



I i i i I i i i i I i i i iTI i 1 1 1 I i i i i I i i i iTi i 1 1 1 I I i i i I i i i i" 



-10 1-10 1-10 1 

log(1+(5 2 ) 



-1 



o 



1 



Figure 10. Environment dependence of the residuals, at fixed M v i r , of the c — Af v j r relation. Different panels correspond to different 
halo mass bins, as indicated (numbers indicate log[M v i r /h — 1 M©]). The solid curve shows the mean, the error bars show the Poisson 
error on the mean, and the dashed curves show the ±lcr scatter. The dashed lines show the overdensity corresponding to a single halo 
in a sphere of radius 2h _1 Mpc with a mass equal to the lower and upper mass limits of each panel. The almost vertical lines have the 
same meaning of Fig. [9] indicating the average overdensity (plus its error) for haloes in a given residual bin. 




Figure 11. Same as Fig. 1101 except that here we show the residuals of the A' — M v i r relation. 
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Figure 12. Same as Fig, 1101 except that here we show the residuals of the q — M v ; r relation. 
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Figure 13. Correlations between the halo variables X' c and c with the galaxy properties So (baryonic disk central surface density, units 
of Mq pc — 2 ), So,* (stellar disk central surface density, units Mq pc~ 2 ), and /^o,fl (stellar disk central surface brightness in the B-band, 
units magn. arcsec -2 ), for haloes with iV v i r > 250, p rms < 0.4, and x a ff < 0.04 (i.e. relaxed haloes). The points are color coded according 
to M v i r (units h -1 Mq) as indicated. 
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Figure 14. Histograms of concentration parameter for different viral masses (units of h~ Mq) and B-band central surface brightness 
(units of magn. arcsec -2 ) for the model galaxies in Fig. 1131 The three mass ranges (from top to bottom) correspond to massive, 
intermediate and dwarf galaxies, while the four surface brightness ranges (from left to right) correspond to low, intermediate, high and 
very high surface brightness. For each subset of haloes the red line shows a gaussian distribution with a mean c and standard deviation 
<j\ n c as given in the top right corner of each panel. Note that the dependence on mass is much stronger than the dependence on surface 
brightness. 



